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■ Abstract. We present various results on the properties of the four infinite sets of the excep- 
tional Xg polynomials discovered recently by Odake and Sasaki [Phys. Lett. B 679 (2009), 
414-417; Phys. Lett. B 684 (2010), 173-176]. These X t polynomials are global solutions 
of second order Fuchsian differential equations with £ + 3 regular singularities and their 
confluent limits. We derive equivalent but much simpler looking forms of the Xi polyno- 
mials. The other subjects discussed in detail are: factorisation of the Fuchsian differential 
operators, shape invariance, the forward and backward shift operations, invariant polynomial 
subspaces under the Fuchsian differential operators, the Gram-Schmidt orthonormalisation 
procedure, three term recurrence relations and the generating functions for the Xg polyno- 
mials. 
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G\ '. 1 Introduction 

O 

Four sets of infinitely many exceptional (Xg) polynomials satisfying second order differential 
equations were introduced recently by two of the present authors [U [2]. They were obtained 
$_i ' as the main part of the eigenfunctions of exactly solvable one-dimensional quantum mechani- 

cal systems which were deformations of the well-known solvable systems of the radial oscillator 
[31 H] and the trigonometric Darboux-Poschl-Teller (DPT) potential [5] by a degree i eigen- 
polynomial. Thus the orthogonality and completeness of the Xg polynomials are automatically 
guaranteed. These polynomials, termed exceptional Laguerre and Jacobi polynomials, have two 
types in each family, LI, L2 and Jl and J2. The Laguerre family LI (L2) is obtained from 
the Jacobi family Jl (J2) by the well-known limit (|E.28j) . which takes the Jacobi polynomials 
to the Laguerre polynomials. The Jl and J2 are mirror images of each other, see (|E.14p . but 
their limits LI and L2 are clearly distinct. These polynomials are exceptional in the sense 
that they start at degree i (£ = 1,2,...) rather than degree constant term. Thus they are 
not constrained by Bochner's theorem [6], which states that the orthogonal polynomials (star- 
ting with degree 0) satisfying a second order differential equations are very limited. Namely, 
they are only the classical orthogonal polynomials, the Hermite, Laguerre, Jacobi and Bessel 
polynomials. 

The concept of exceptional orthogonal polynomials was introduced in 2008 by Gomez-Ullate 
et al. |H]. Within the Sturm-Liouville theory they constructed X\ Laguerre and X\ Ja- 
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cobi polynomials, which turned out to be the first members of the infinite families. The re- 
sults in [H] were reformulated in the framework of quantum mechanics and shape-invariant 
potentials [9] by Quesne and collaborators [101 111) . They found the first member of the de- 
formed hyperbolic DPT potential family, which was also given in [I]. Quantum mechanical 
reformulation offers two merits. Firstly, the orthogonality and completeness of the obtained 
eigenfunctions are guaranteed. Secondly, the well established solution mechanism of shape in- 
variance combined with Crum's method |12j . or the so-called factorisation method [3J or the 
susy quantum mechanics [4] is available. For systems consisting of discrete eigenvalues only, 
shape invariance is a well-known sufficient condition for exact solvability of one-dimensional 
Schrodinger equation. The discovery of the four sets of infinitely many exceptional orthogonal 
polynomials was achieved by pursuing shape invariant deformation [HE]. To be more pre- 
cise, deformation of the potential by keeping the property of shape invariance. After the first 
paper on infinitely many Xn polynomials [T], Quesne reported a type II X<i Laguerre polyno- 
mials [131 114) . which led to the discovery of the L2 family of Xg Laguerre polynomials [2]. In 
a previous paper [15], two of the present authors unveiled infinitely many polynomial identi- 
ties of degree 31 involving cubic products of the Laguerre or the Jacobi polynomials, which 
encode the information of exact solvability of the differential equations governing the Xg poly- 
nomials. 

In this paper we explore various properties of the Xi polynomials. We emphasise that the Jl 
and J2 polynomials are the global solutions of a Fuchsian differential equation having i+3 regular 
singularities. They are located at ±l,oo and the i zeros of the polynomial &(t7; A) ()2.2j) . (|2.4j) 
which is used for the deformation. Factorisation and shape invariance are reformulated ac- 
cordingly, leading to Rodrigues formulas and the forward and backward shift operations. The 
existence of the extra regular singularities implies that the ordinary vector space spanned by 
degree n polynomials Vo, n = Span [1, x, . . . , x n ] is not invariant under the Fuchsian differential 
operator (|6.2p . Appropriate invariant polynomial subspaces are introduced and their properties 
are used to derive the explicit forms of the exceptional polynomials. Some of the important 
subjects in orthogonal polynomial theory, namely, the Gram-Schmidt orthonormalisation, the 
generating functions, three term recurrence relations, the zeros of these orthogonal polynomials, 
etc, are also discussed. New infinitely many polynomial identities underlying the forward and 
backward shift operations are also reported. 

The plan of this paper is as follows. In Section[2]we present the explicit forms of the four sets of 
infinitely many exceptional orthogonal polynomials together with their weight functions and the 
normalisation constants. They are equal to those reported earlier [H[2], but look much simpler 
than the original ones. The new forms of the polynomials reveal the structure of the theory. In 
Section [3] the Fuchsian differential equations governing these polynomials are discussed. Shape 
invariance and Rodrigues formulas are presented in Section The identities underlying the 
forward and backward shift relations are presented in Section [5] The polynomial subspaces 
invariant under the Fuchsian differential operator are discussed in Section [6) This gives another 
concise proof of the new forms of the exceptional polynomials. In Section [7] we provide the 
integration formula which is essential for the Gram-Schmidt construction in Section [HI Section [9] 
gives the generating functions for the Xg polynomials. The double generating function, that is, 
the generating function of the generating functions, is presented for the LI and L2 exceptional 
Laguerre polynomials. A substitute of the three term recurrence relations for the Xg polynomials 
is introduced in Section [TUJ In Section [11] we state the qualitative features of the extra zeros 
of the Xn polynomials without proof. The final section is for a summary and comments. Some 
technical details are relegated to Appendices. The equality of the new and original forms of 
the Xn polynomials is demonstrated in Appendix [A] Forward and backward shift relations are 
proved in Appendix [B] Derivation of the integration formula is provided in Appendix [C] The 
properties of the Xg Jacobi polynomials as the solutions of the quantum mechanical systems 
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with deformed hyperbolic DPT potentials are summarised in Appendix [Dj A concise summary 
of some important properties of the Laguerre and Jacobi polynomials is given in Appendix [E] 
for self-containedness. 

Throughout this paper we stick to the notation of our previous papers [H [21 [151 [16] . Reference 
to the quantum mechanical language is made minimal in order to attain wider readership than 
before. Most concepts and formulas are common to the four sets of exceptional polynomials. As 
far as possible we use generic formulas valid for all the four different sets of Xi polynomials, in 
order to emphasise the underlying structure and at the same time to avoid redundancy. 



2 Exceptional Laguerre and Jacobi polynomials 

Here we present four sets of infinitely many exceptional orthogonal polynomials [H [2], among 
them two are deformations of the Laguerre polynomials, and the others are deformations of the 
Jacobi polynomials. They are expressed as a bilinear form of the original polynomials, the 
Laguerre or Jacobi polynomials and the deforming polynomials, depending on the set of parame- 
ters A and their shifts S and a non-negative integer £, which is the degree of the deforming 
polynomials. It is important to stress that these explicit forms (|2.ip - (|2.4p can be derived by 
the Darboux-Crum transformations starting from the original (Laguerre or Jacobi) polyno- 



mials \17\ 118]. The two sets of exceptional Laguerre polynomials [l = 0, 1, 2, 



n 



0,1,2,...) 



arc: 



A + 8)P n {i T , g + £-!)- &( m X)d v P n ( V ; g + £ - 1) 
(n + g + \Y X {{g + \)Um * + S)P n (r); g + £ + 1) 



+ r J ^(r J ;\)d r] P n (r J ;g + £ + l)), 



(LI) 
(L2) 



(2.1) 



in which A = f g > and 5 = f 1 and 



dcf I 



(-V), (LI) 



4 9 ' (L2) 



(2.2) 



The two sets of exceptional Jacobi polynomials (£ = 0,1,2, .. ., n = 0, 1, 2, . . .) ara^: 



' ( n + h + 1) 1 ((/> + I)^; A + S)P n (rj; g+ £-l,h + £ + l) 

+ (1 + VMV, X)d v P n (v, g + £-l,h + £ + l)), (Jl) 
n + g + \)~\(g+ \)Um A + S)P n ( mg + £ + l,h + £-l) 

- (1 - 77)6(77; X)d 11 P n (r ] ; g + e + l,h + £-l)), (32) 



(2.3) 



in which A = f (g, h), g > 0, h > 0, S = f (1, 1) and 



P n {mg,h)= Pt hh ~ l \i), 



P ( e 9+e 1 h 1 *\t,), g>h>0, (Jl) 

(-g-t-kW-*)^ h>g> ^ (J 2 ) 



(2.4) 



It should be remarked that the naming Jl and J2 are interchanged from the previous ones [2] in order to 
respect the logical consistency rather than hysteresis. 
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It should be stressed that the deforming polynomial &(rj; A) does not have a zero in the domain 
of orthogonality, (0, oo) for LI and L2 and (—1,1) for Jl and J2. This is demonstrated explicitly 
in (2.39), (2.40) of [15]. These Xg polynomials have the following general structure 



do(n, \)P e , n (m X) = di(A)&(»7; A + S)P n (r]; \ + iS + S) 

- d 2 (vMv, X)dr,Pn{m x + es + s), 



(2.5) 



where 



d (n, A) 



clef 



n + g + \, 
n + h + \, 



di(A) 



def 



1 
2' 



2' 



(LI) 
(L2, J2) 

(Jl) 

(LI) 
(L2, J2) 

(Jl) 



rf 2 (r/) 



def 



def J + 1, 


(L1/L2) 


= \t(i,-i), 


(J1/J2) 


'1, 


(LI) 


< -7], 


(L2) 


,T(l±r?). 


(J1/J2) 



(2.6) 



We introduce a linear map for a differentiable function p(r]), 

[?(»/)] = f di(X)&(m X + <5)p(tj) - d 2 (7])£ e (? r , X)d r ,p{ri). 
Then the A^ polynomial (|2.5[) is expressed succinctly as 
do(n, AJP^fa; A) = E <)A [P n fa; A + M + 5)] . 



(2.7) 



(2.c- 



The Xq polynomials Pq n (r/; A) = P n (r]] A) are the undeformed polynomials, i.e., the Laguerre 
or the Jacobi polynomials themselves. Therefore the above formulas (|2. 1|) . (|2.3[) and (|2.5p for 
£ = are non-trivial identities among the Laguerre or the Jacobi polynomials 

d (n, X)P n ( V ; A) = di(A)P n fa; A + 6) - d 2 (r ] )d 71 P n (i 1 ; X + 8). 

This is shown by using (1R21) . (IETToT) and the following; (lETTTjl for LI, (|ETTT|) and (IETT21) for L2, 
(IE.22P for Jl and (IE, 25ft for J2. The Xi polynomials Pi, n {j}\ A) are degree l + n polynomials in rj 
and start at degree I: 



Pto(v,X) = &(ri;X + 5). 



(2.9) 



They are orthogonal with respect to the weight function V\^(r?, A) which is a deformation of the 
weight function W(rj; A) for the Laguerre or Jacobi polynomials: 



/ 

m(v,x) 



Pl,n{rr, X)P£, m (v, X)Wi(r]; X)dr] = h^ n {X)5 n 
def W(ri; X + IS) 



(2.10) 



W( V ;X) 



def 



Uv,x) 2 



< i] < oo, 

-1 < 7] < 1. 



(L) 
(J) 



, ^(1-^-5(1 + ^3, 
The normalisation constants are meromorphic functions of the parameters g, h and 

def 1 



(2.11) 



(L): K(g) 



2nl 



T(n + g + ±), 



(2.12) 
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( n+g+2£-l , . 



htM^hnb+e)* (2.13) 

(J) : nn + 9+ l)T(n + h + ^ 

y ' vy ' 2n\{2n + g + h)Y{n + g + hy y ' 



he, n (g, h) = f h n (g + £,h + £) x < 



f (n+fe+l+i)(n+g+2l-±) 
(n+/i+i)(n+g+f-i) ' 

(n+g+l+\)(n+h+ 2 l-\) 
(n+ 9 +|)(n+/i+£--|) " 



(Jl) 

(J2) 



(2.15) 



As stressed in § 3 of [2], the Jl and J2 polynomials are the mirror images of each other, in 
the sense r] —77 and g -H- h, as exemplified by the relation £$ 2 (r);g,h) = (—iy£f l (—r);h,g). 
However, they lead to the two different sets of the exceptional Laguerre polynomials, Jl— )-Ll, 
J2— >L2. In terms of the limit formulas 

hmP^(l-f)=4*)(±x), (2.16) 

p— too ^ 

it is easy to see the relations P^ n (ry;A) (|2.3p — 7- (12.ip together with the deforming polynomials 
&(ry;A) (|2^]l -)' ([2^p . the normalisation constants [ (12~TT?D with ([235]> ]-»[ ([210]l with ([27T3]) ] and 
others. The explicit forms of the polynomials f)2. 1 1) and (j2.3[) are much simpler than those 
given in the previous papers [HE]. In Appendix [X] we will give simple demonstration that these 
apparently different forms of polynomials are in fact equal. 

The basic ingredients of the theory of exceptional orthogonal polynomials are the base poly- 
nomial P n and the deforming polynomial ^. They satisfy the second order differential equa- 
tions (IE.4P and (IE. 17ft . which can be expressed as 

csjfajagPnfa; A) + Cl (?y, X)d n P n { m X) = -\S n {X)P n { m A), (2.17) 
c 2 (r l )^ i (v,X) + cx(v,\i)d^i(v,><) = -\£e(XM V ;X), (2.18) 

where 

U ' \h-g-(g + h + l) V , (J) jl- ? y 2 , (J) 

c^A,W ±( / + '~ l + ??) ' 

\ T ( g + h + 2£-l + (g-h)ri), 

£ n (X)^{ 4n > (L) ^(A)^J^ ( 2 .i 9) 

V ' \4n(n + 5 + /i), (J) K ' \4£(£±gTh-l). (J1/J2) 

The deforming polynomial ^(ry; A + 5) is expressed in terms of &(ry; A), 

di(A)&fa; X + S) = d 1 (X + £8)U m X) + d 2 {v)d n Um A), (2.20) 

with <ii and d 2 defined in (|2.6p . This is shown by (|E.2p . ()E.15p and various identities of the 
polynomials; (jETTTj) for LI, (jETTTj) and (lETT2|) for L2, (|ET22l) for Jl and (|R25l) for J2. Conversely 
^(ry; A) is expressed in terms of &(ry; A + 5), 

d 3 (A + 4$, ^(ry; A) = d 3 (A, £)Uv, X + 5) + ^\ d,Uv, A + S), (2.21) 
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where 



d 3 (X,i) 



def 



g + t 
1, 

h + £ 



i 

2> 



1 

2 ' 



(LI, Jl) 

(L2) 

(J2) 



(2.22) 



This is shown in similar ways as above. The Laguerre and Jacobi differential equations for P n 
(|2.17p can be factorised into the forward and backward shift relations for P n : 



c T d v P n (rj; A) = / n (A)P n _i(r?; A + 6), 

ci(v, X)Pn-i{v\ A + 5) + c 2 (n)d v P n _ 1 (n; A + 6) = -\cj,b n ^i{X)P n {n; A) 



where 



def I 2, 



(L) 
(J) 



/n(A) 



-2, 

-2{n + g + h), 



(L) 
(J) 



6„_i(A) 



-In. 



(2.23) 
(2.24) 



(2.25) 



See (|E.2p . (|E.3[) and ()E.15p . ()E.16p for the explicit forms of the forward and backward shift 
relations. Note that f n (X) and 6 n _i(A) are the factors of the eigenvalue 



£nW = /n(A)6 n -l(A), 



n = 0,l,... 



3 Fuchsian differential equations 
with extra £ regular singularities 

The exceptional Laguerre and Jacobi polynomials satisfy a second order linear differential equa- 
tion in the entire complex rj plane: 

n e (x)p itn ( V ; a) = s t , n (\)Pt, n (m a), s t , n {\) = £ n {x + ts), (3.1) 

in which the eigenvalue £ n is defined in (|2.19p . The polynomials are not constrained by 
Bochner's theorem [B] by the very fact that they start at degree t (|2.9p instead of degree 
constant term. As with the Laguerre and Jacobi differential equations (j2. 17j) . the second order 
differential operator He(X) is factorised into the product of the forward shift operator J^(A) 
and the backward shift operator Bg(X): 

H e (X)= B e (X)Ti(X), (3.2) 

, t(x) ^ Cr i^^_ dMX+S) y (3 . 3) 



The differential operator T~Le(X) is obtained from the factorised quantum mechanical Hamilto 

02 fa) 



nian (16) of [T] and (3) of [2], by similarity transformation. Note that EliZZ^+jj) = \ g W(n; A+ 



+ 1)S). It is straightforward to derive the explicit form of T-Li{X): 
He(X) = -4 (c2(ri)^ + ( Cl (r/, A + 18) - 2c 2 (n)d v log A)) 

+ Ml(A) « + 4 + 5) • (3 - 5) 
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Use is made of the second order differential equations for £^(77; A + 6) ()2. 18[) and the identi- 
ty (|2.20p to derive the above simple result. 

For £ = the above differential equation (|3.ip with (|3.5p reduces to the second order 
differential equation for the Laguerre or Jacobi polynomials: 

«o(A) = -4 (c 2 ( V )^ + Cl (r?, X)-f] , (3.6) 



<i?7 2 ' dij 

(L): ^(A) = -4^ + ( 5 +|-r ? )A), (3.7) 

vd%Pn(v,g) + {a + \ -v)9 v Pn(rr,9) + nP n {ma) = o, 

(J): £ (A) = -4 + s-fc + fc + i^iLy (3.8) 

(1 - V 2 )d*P n (m g,h) + (h-g-(g + h + 1)77)^^(7?; /i) 
+ n(n + 3 + h)P n (r)\ g, h) = 0, 

which has, as is well-known, one regular singularity at 77 = and one irregular singularity at 
77 = 00 for the Laguerre case and three regular singularities at 77 = ±1, 00 for the Jacobi case. For 
a non-negative integer £, the singularity structure of the second order differential equation ()3.1|) 
is again quite simple. It has extra regular singularities at the £ zeros of the deforming polyno- 
mial ^(77; A): 

? ? = 7 ?i> 6(%;A) = o, j = 1,2, ...,£, 

and the corresponding exponents are the same for all the singular points: 

{exponents at rjj} = {0,3}, j = 1, 2, . . . , £. 

In other words, the alternative solution of the second order linear differential equation (|3.ip has 
a cubic zero at 77 = rjj. It is singular at r/ = 0, 00 for the LI and L2 and at 77 = ±1 for the Jl 
and J2 with the exponents replaced by g — > g + £, h — > h + £. To the best of our knowledge, 
the two sets of Xi orthogonal polynomials, Jl and J2 are the first examples of global solutions 
of Fuchsian differential equations having as many as £ + 3 regular singularities and forming 
a complete orthogonal system. The LI and L2 are confluent types obtained from Jl and J2 
by certain limits (IE.28p . By the way, it is elementary to show that £0(77; A) has only simple 
zeros. 

Let us write down the explicit form of the above differential equation for the four cases: 
(LI) : V d*P e , n (r,; A) + (g + £ + \ - 77 - 2 ^f|^) Wfa A) 

+ (2 UmX) +n-£j PUV, A) = 0, (3.9) 
(L2) : T$P t)n ( m A)+ (^g + £ + \- r) - 2^|^) d.P^r,; A) 
1 ( (g + ^£1(77; A + <5) \ 

+ 1-2 ^(77- A) h n + ^ I Pe^niv'i A) = 0, (3.10) 

(Jl): (l-^P^faA) 

+ (fc - 5 - (<7 + * + 2i + 1)77 - 2 (1 "^f ;X) ) VU* A) 
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/ 2(h + ±)(l-r))d v &( V ;\ + S) , h , f\\ 
+ t-, — rr h £{£ + g - h - 1) + n(n + g + h + 2£) 

\ / 

xP<, n (i);A)=0, (3.11) 
(J2) : {l- V 2 )d 2 v P e , n ( m \) 

+ (h-g-(g + h + 2i + l)r,- 2 (1 "^f ;A) ) <VWtt A) 

+ — \- £(£ + h — g - 1) + n(n + g + h + 21) 

\ J 

xP e>n ( m X) = 0. (3.12) 

Let us remark that the zeros of the shifted deforming polynomial £^(77; A + d) are regular points. 
It is straightforward to verify by direct calculation that the LI, L2, Jl and J2 polynomials 
(|2.ip - (|2.4p for lower £ and n really satisfy the above differential equations (|3.9p ~ (|3.12p . For 
analytical proof see the subsequent sections. 



4 Shape invariance 

This section is a reformulation of the shape invariance in the language of ordinary differential 
equations. In one-dimensional quantum mechanics, shape invariance is a sufficient condition [9] 
for exact solvability and it was the guiding principle for the discovery of these Xg orthogonal 
polynomials [U [T6| [T5| E|. Let us introduce another second order linear differential operator 
(A) by interchanging the order of the two factors Fg{X) and Bg{X) of %e{X) (|3.2p : 

n1\\)= T,(X)B £ (X). 
It is obvious that these two operators are intertwined by ^(A) and Be(X): 

Fz(X)B e (X)Tt(X) = Te{X)Ut{X) = H?\x)T e (X), 
B e (X)T e (X)B e (X) = H e (X)B e (X) = B e (X)nf ] (X), 

which implies that these two associated linear differential operators T-Li(X) and (A) are 
iso-spectral except for the lowest eigenfunction P^ 0(77; A) (|2.9p which is annihilated by J-^(A), 
see (pj3j) : 

TtWttiv, X + S) = (=> Ht(\)P t ,o(m A) = 0) . 

If we denote the set of eigenfunctions of (A) as {P^„ '(77; A)} {n = 0, 1, . . .), with arbitrary 
normalisation, we obtain one to one correspondence of {Pg in {v'i A)} and {P^(rj; X)} except for 
the lowest eigenfunction of Hg(X): 

He(X)Pt,n(iK A) = £i,n{X)Pi, n (m A), n = 0, 1, . . . , 

n^WP^imX) = s tjn (X)p$_ 1 ( m x), n = 1,2,..., 

^(A)?^^; A) cx P^-ifa A), B£(X)P^_ 1 ( r n; A) cx P^fa; A). (4.1) 

This much is a trivial consequence of the factorisation of He(X) f|3 . 2 [) . The essential property of 
the two associated differential operators Hi(X) and 'H < f\x) is that T~6p {X) has the same shape 
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as 7{e(X) with shifted parameters, A — > A + 8, and an additive constant corresponding to the 
lowest eigenvalue measured from the bottom: 

Hf\\) = H i (X + S)+£ 1 (X + £S), or 

T t {X)B t {X) = B £ (X + 8)T t {X + S) + £ 1 (X + £8). (4.2) 

The system is called shape invariant if the above condition is satisfied. See [H [2j [15] for the 
definition of shape invariance within the framework of quantum mechanics. The above is an 
equivalent definition within the framework of the Sturm-Liouville theory. 

Shape invariance allows us to choose the normalisation of {P^(r/; A)} to achieve 

fa; A) = P e>n ( m A + 8), n = 0, 1, . . . . 

Then the above two relations (|4.ip . the forward and backward shift relations, give the constraints 
on the functional forms of {P^ n (r/; A)}: 

T e (X)P e>n (r ] ;X) = f n (X + £8)P e , n . 1 ( V ;X + 8), n = 0, 1, . . . , (4.3) 
B i (X)P itn - 1 ( m X + 8) = b n - 1 (X + e8)P itn {r 1 -X), n = l,2,..., (4.4) 

where f n (X) and 6 n _i(A) are given in (|2.25p . These amount to a version of Rodrigues formula 
expressing Pg tTl (ri;\) in terms of repeated application of the backward shift operators on the 
lowest degree eigenpolynomials with the n-th shifted parameters, P^ofa; A + nS) = A + 
(n+l)8): 

k=o (£ + k)S) 

n-1 

where Y\ a k = clqcli • • • a n -i. F° r the Xg Laguerre and Jacobi polynomials, the Rodrigues 

k=0 

formula reads explicitly: 

L : P t n (i ; A = — — — M ( -j- + log — — 

n! e-Vrf+t-2 \drj ^(r/; A + fed) / 



&(rr,X + n8) 



(i\ P ( \\ (-1)" frfaA) "rjfd ^ ^(r,;X + (k + l)S) \ 

(J): ^M) = W(i _ r4(i + r4 n^ + ^lo g UvX + kS) ) 



77)^(1+^-^=0 
} ^ : ( ' + ;; -6 fa A + (n + l)d). 

For i = 0, ^ofai A) = 1, the above two formulas reduce to the well-known Rodrigues formulas 
for the Laguerre and Jacobi polynomials, (|E.6|) and (|E.19p . For £ = 1 these two formulas are 
equivalent to the Rodrigues-type formulas (77) and (52) in Gomez-Ullate et al. [7j. For lower £, 
it is straightforward to verify the shape invariance relation (|4.2p by direct calculation. In [15] 
it was shown that the shape invariance relation is attributed to a new polynomial identity of 
degree 3£ involving cubic products of the Laguerre or Jacobi polynomials. These identities are 
proved elementarily by combining simple identities in [15j . 

It is a good challenge to derive more explicit expressions of the Rodrigues formulas for each 
of the four families of the exceptional orthogonal polynomials. 
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5 Forward and backward shift relations 

Again it is straightforward to verify that the above Rodrigues formulas provide the explicit forms 
of the four families of the Xg polynomials (|2.ip ~ (|2.4p for lower £ and n by direct calculation. For 
analytic proof, one needs to verify that the above forward (|4.3|) and backward (|4.4|) shift relations 
are actually satisfied by the four families of the Xg polynomials (|2.ip - (|2.4p . The forward and 
backward shift relations are attributed to new polynomial identities of degree 2£ + n — 1 and 
21 + n involving cubic products of the Laguerre or Jacobi polynomials in a similar way to the 
shape invariance relations |15] . The forward (F) (14. 3p and backward (B) (|4.4p shift relations 
read explicitly 



(F) : = c T ( A + 8)— - d^eim X + S)) Pg^rj; X) 



d 

dr/ 

- f n (X + l8)Um X)Pi, n -x(w X + 5), (5.1) 
(B) : 0= (^ e (r,;X)(c 2 ( V )-^ + c 1 (r l ,X + £d)) - c 2 (r l )d r ,&( m X)\ P^ir)-, X + 6) 

+ \c T b n ^{X + £5)Um x + S)P t , n (m (5-2) 



where P^ n is given in (|2.5p . 

The new polynomial identities are as follows: 
(LI): ( a = g + £-\) 

(F) : = (Lf\-x)l - d xL f\-x)) {Lf\-x)L^\x) - L^H-xfeL^H*)) 

+ 4 a ~ 1] {-x) (x) - (-ziO.L&ix)) , (5.3) 

(B) : = (Lf~ l \-x)(x± + a + 1 - x) - xd x Lf- l \-x)) 
x (L^i-x^Six) - Lft-!)^)) 

- nLf ] (-x) {4 a) (-x)L^ (x) - {-x)d x L<T* (x)) , (5.4) 
(L2): ( a = g + £-\) 

(F): 0=(L ( f a - 2 \x)±-d x L ( f a - 2 \x)) 

x {{a-£ + l)L { - a - 2 \x)L^ +1 \x) + xL { ~ a ~ 1 \x)d x L i ^ +l \x)) (5.5) 
+ 4" Q - 1 )(x)((a - £ + 2)L ( f a ^\x)L^ 2 \x) + xL^ a - 2 \x)d x L^ 2 \x)) , 
(B) : 0=(L ( f a - 1 \x)(x£+a+l-x)-xd x L { - a - 1) (x)) 

x((a-£ + 2)L ( f a - 3 \x)Lt\ 2) (x) + xL^ 2 \x)d x L^\x)) (5.6) 

- nL^ a ~ 2) (x) ((a - £ + 1)4" Q " 2) {x)L[^ (x) + xL { f a ~ 1] (x)d x L^ (x)) , 

(J2): {a = g + £-\, p = h + l-\) 

(F) : 0= (pj-^^-^Pj^^^ 

- (1 -x)Pt a - 1 ^- 1) (x)d x Pt +1 ' (3 - 1 \x)) - fr + a + P + l)Pt a ~ 1,P ~ 1) (x) 
x((a-£ + 2)Pt a -^ +l) (x)P^ + 1 m (x) 

- (1 - x)Pt a - m {x)d x Pt\^\x)), (5.7) 
(B) : 0= (p^ a ~ 1 ^ 1 \x)((l-x 2 )£ + P-a-(a + (3 + 2)x) 
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- (1 - x 2 )^]-*-^- 1 ^)) ((a-£ + 2)P^ +l \x)Pt + 1 m {x) 

- (1 - x)P^-^\x)d x Pt\ 2 '"\x)) + 2nPt a - m (x) 
x((a-£ + l)Pt a ~ 2 ^ ) (x)Pt +1 ^- 1) (x) 

- (1 - x)P { - a -^- 1) (x)d x Pi a+ ^- 1 \x)). (5.8) 

We do not show those for the Jl polynomials, since Jl and J 2 are related by the mirror image, 
rj -H- — T] and renaming of the coupling constants g <->• h. Thus the corresponding identities are 
essentially the same. The derivative terms in the above identities (I5.3p -( j5"78j) can be eliminated 
by using the forward shift relations (|E.2p and (|E.15P for the Laguerre and Jacobi polynomials. 

In Appendix El the forward and backward shift relations for the four types of Xi polynomials 
are proven elementarily. 



6 Invariant polynomial subspace 

Another characterisation of polynomials satisfying differential equations is the existence of in- 
variant polynomial subspaces |8j. For the Laguerre and Jacobi differential equations with T-Lq(X) 
(]3.7p . (I3.8p . the space Vo n of degree n polynomials in r\ is invariant 



■Ho(A)V , n Q Vo, n , Vo, n = f Span[l,7/, . . . ,rj n ], 

fi (X)r] n = £. n {\)ri n + lower orders. (6.1) 
However, Vo, n is obviously not invariant under Tlii{X) (|3.ip . (|3.9p - (|3.12p : 

W*(A)V ,„ 2 V 0>n , £ = 1,2,.... (6.2) 
Instead we have 

HtWV t , n C V<,„, 1 = 1,2,..., 

def (Sp&n[r j k &{r] ] g + l)-kT ] k ~ 1 &(r ] -,g);k = 0,l,...,n\, (LI) 
Vln = \ (6.3) 

[Span[r, k ((g + ±)&( ?? ; g + 1) + fc&fa; 5 )) ; fe = 0, 1, . . . , n] , (L2) 



Span[? ? fc (/i + h)Zt(w, 9 + l,h + 1) 

+ k(l + r/)^- 1 ^^; g, h); k = 0, 1, . . . , n] , (Jl) 
Span[n k (g + ^ e (r];g + I, h + 1) 

- k(l - rirf-^fag, h); k = 0, 1, . . . ,n] . (J2) 



(6.4) 



As the basis vectors of the invariant polynomial subspace Vn n for the Xi Jacobi polynomials, 
one could have chosen for k = 0,1,... , n, 

(l-7 ? ) fc - 1 ((/ i +i)(l-7 ? )^(r ? ; 5 + l,/i+l)-^(l + ^(^5,/i)), (Jl) (6.5) 
(l + Tj) fc ((fc+±)&fa; g + l,h + l) + k £ t {r};g,h)), (Jl) (6.6) 

(l-r)) k ((g + lM mg + l,h + l) + k&(r)-,g,h)). (J2) (6.7) 

After overall rescaling, (I6.5P and (|6.7[) go to those of the Xi Laguerre polynomials (|6.3p in the 
limit (|E.28p . as Jl— >L1 and J2— )-L2. As will be shown shortly these basis vectors (|6.6p and (j6.7j) 
have simpler integration formulas than (I6.4H . It should be stressed that these basis vectors have 
the same common structure as the Xg polynomials (|2.8p . ([27 



Sf,A[pfc(f?)] = <MA)&(??;A + $)Pk{v)- d 2 {r])it{r];X)d ri p k (i]) G V^ n , fe = 0, 1, . . . , n, (6.8) 
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in which pk(i]) is an arbitrary degree k polynomial in rj. These basis vectors are so chosen as not 
to develop any singularities at the zeros of £e(rj; A) when applied by the forward shift operator 
Fz{\) (|3.3p . In fact, for a polynomial p(rj), Te(X) acts on Et t x\p(v)] as 

c-^tiX^xlpirj)] = {d 1 (X + 8)d v p( V )-d v (d 2 ( V )d v p( V )))( e ( V ;X + 8) 

+ d 2 (r j )d r ,p(r ) )d ri ^(r ) ;X + S), (6.9) 

where we have used (|2.20p to eliminate 9„^(ry; A). Since the backward shift operator Bi(X) (j3.4|) 
does not cause any singularity at the zeros of Ciiv'i ^) an d the operator T-Ln{X) = Bi(X)J r e(X) has 
the form ()3.5p . the application of T-Le(X) on these basis vectors will result in polynomials in rj. 

Let us evaluate the action of the second order differential operator He(X) on E£ t x[p(ri)] by 
applying the backward shift operator Bg(X) on (I6.9P : 

He(X)E ij Mv)] = Uv, A)X(r?) + &(r/; A + S)Y( V ). 

The coefficients X(r]) and Y{rf) can be written as 

X( V ) = 4d 2 (ri)d v (c 2 (r ] )d 2 v p(ri)+c 1 ( V ,X')d v p(r ] )) = -d 2 ( V )d v {n (X')p( V )) , 
Y(rj) = -4di(X)(c 2 (ri)d^p(i]) + c\(i], X')d v p(rj)) = di(X)H (X')p( V ) , 

where A' = A + £8 + 8 and Ho(X) is defined in (|3.6p . Therefore the results can be expressed in 
quite a simple form as 

H e (X)E ejX \p(r))] = ~ e>x [H (X + £8 + 8)p( v )] . 

By taking p(rj) = j>fc(ry), the above basis vectors satisfy 

Hz(X)Ee t xMv)} = ^x[H (X + 18 + 8)p k ( V )] e V e , n , k < n. 

In other words, we have shown that V^ n is an invariant polynomial subspace of the differential 
operator He(X) and we obtain, corresponding to (16. If) . 

Ht(X)E£ t x\p n (.v)] = ^n(A + £8)Ee : x\Pn(v)} + low er orders. 
In particular, if we choose p n {r]) as the eigenfunction P n (r/; X + £8 + 8) of Ho(X + £8 + 8), 

■H (X + £8 + 8)P n (rr, X + £8 + 8) = £ n {X + £S)P n ( V ; X + £8 + 8), 
then we find that E^x[Pn(Vj X + £8 + 8)] is the eigenfunction of "H^(A), 

H e (X)E e> x [P n {m X + £8 + 8)]= £ n (X + £8)E e , x [P n (m X + £8 + 8)], 

as given in (|2.8p . Note that the effect of 8 in the eigenvalue cancels out. This is another 
analytical proof for the explicit forms of the polynomials P£,n(Vi ^) as given (|2.8p . 



7 Integration formulas 

Another well-known construction method of orthogonal polynomials is the Gram-Schmidt or- 
thonormalisation of certain basis vectors under a given inner product specified by a weight 
function. Let us introduce two types of inner products (*, *)e t \ and (*, *)x- 

(p(r]),q( 7 l))^x = / Pirfqi^WeivA)^, {p(v),q{v))x d = / P(v)q(v)W(rj;X)dr], 
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in which p{rf) and q{rf) are arbitrary functions and We(r]; A) is the weight function for the 
Xi polynomials, whereas W(rj; A) is the weight function for the Laguerre or Jacobi polyno- 
mials (f2TL|) . 

Here we present the integration formulas: 

= d l (X)d 3 (X + £d,e)(p(r l ),q(r l )) x , + l4(d vP ( V ),d v q( V )) y+d , or (7.1) 
Z £t x[p{v)] Hf,A[?('7)]>Vi(77; \)dr) 

di(A)d 3 (A + IS, £)p(v)q(v)W(v, V) + \c% d^d^Wirj; A' + 8)y v , (7.2) 

where A' = A + £8 + 8. For a proof, see Appendix O 

With this formula one can easily verify the orthogonality relation and the normalisation 
constants of the X^ polynomials (|2.10p . (|2.13p . f|2. 15|) given in Section [2j One simply takes 
p(rj) = P n (T}', A') and q(rj) = P m (rj; A') for L1-J2. Then the two terms in (|7.2p read 

(P n (r ] ;\'),P m ( m X')) xf = J P n (r ] -X')P m (r 1 -X')W{r l -\')dri = h n (\')5 nm , 

4 J a v p n { m \')d v p m ( m \')w(m x ' + s)dq = iW(A')X-i(A' + 8)s nm , 

where we have used (|2.23[) . Therefore, from (|2.8p . the normalisation constant hi n is expressed 
as 

d (n, \) 2 h^ n (\) = di(A)d 3 (A + £8,£)h n (\') + ±/ n (A') 2 /i„-i(A' + 8), 

where A' = A + 18 + 8 and h-\{\) = f 0. This formula with h n (12J2]) and (gUj) gives (pT3l) 
and (^T5|) . 

The integration formula reveals another important property of the invariant subspaces, the 
orthogonality, Ve, m 1 P£,„(r?; A), m < re: 

(S^ jA [p m (77)], P^„(r?; A))^ A = 0, m < n, 

which is a simple consequence of (|2.8p and the well known fact 

(p m (?7),P„(r?; A)) A = 0, m < n. 

8 Gram— Schmidt orthonormalisation 

As for the direct application of Gram-Schmidt orthonormalisation, one orthonormalises 
{z,£,\\pn(v)]} i n Q6.8j) with respect to the inner product (*, *)^,a- The following choice of the 
function p n is made: 

Pn(v) = I 77 "' (8-1) 

Pn{ " \(i±r?r, (ji/j2) 1 ' 

with which the inner products (p n (v) > Pm{ r i)) \' axe easily expressed in terms of gamma functions. 
This corresponds to the choice of the basis vectors in (|6.6p and (|6.7p . 
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The Gram-Schmidt orthonormalisation procedure at the n-th step determines P^ jn (r/; A) by 
the following formula 

n-1 

Pt >n (ri;X) oc E i:X [pn(v)] ~ ^2 P t,miv, A)<P^, m (??; A), Ei yX [p n {ri))) ^ x h^ m (X)~ l , 

m=0 

which is to be compared with the procedure for the undeformed {I = 0) polynomial 

n-1 

P n {rr,\') CXPniv) ~ Yl ^'){Pm{V\ ^')iPn(??)) A /^m(A')~ 1 , (8.2) 

m=0 

where A' = A + 18 + 8. These two orthonormalisations are essentially the same. In fact, by 
applying E^[-] to (j8.2|) and using (|2.8p . we have 

n-1 

Pt,n{lT, A) OC Ei t x\p n (v)] - Y P t,m(V> X)d*o(m, *){Pm{m A / ),Pn(??)) A //lm(A / ) _1 . 

Comparing these we obtain 

= d (m,X)(P m (r];X'),p n (r])) yh^X'y 1 , m = 0,...,n- 1. 
In other words, the Gram-Schmidt orthonormalisation for the undeformed (I = 0) polynomials 

Pn(v) -> -Pn(??; A'), 
provides that for the Xi polynomials, too. 



9 Generating functions 

Generating functions for orthogonal polynomials have played another important role in classical 
analysis. Let us define generating functions for the polynomials Pg^ n and for the undeformed 
polynomials P n , 

oo oo 

G e (t, r,; X) d ^ f t n Pi,n(v, A), G(t, V ; A) t f G (t, 77; A) = £ t n P n ( m A). 

n=0 n=0 

The latter is quite well known 



G(t, V ;X) 



'G^(t, V ), (L) 
G^h h -^(t,r]), (J) 



where G^Xt^x) and G^(t, x) are given in (|E.8[) and (|E.2ip . Since P^ in (r/;A) is expressed 
linearly in terms of P n (r/; A') (|2 . 5[) (A' = A + £<5 + 8), the generating function Gi(t,r];X) is 
expressed simply in terms of the known ry; A'): 

00 

do(*$, A)G,(t, ??; A) = tn Mn, \)Pi,n{m A) 

n=0 

= di(A)&fa; A + 8)G(t, 77; X') - d 2 (r,Mv, \)d v G(t, m X'). 
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The forward shift relation (|2,23p implies 

d v G(t, rj; A) = c-H f tdt+1 (X)G(t, m A + 5). 
Here we present the concrete forms of the generating functions: 



do(tdt,\)Gi(t, m \) 



IL 



-V) + Th L 



{(9 + h)L 



{r}) + 



(1+V)t t 1 + 9+f_ 



R \R 1 1+R-t 1 1+R+t 



+ 



{h+\)P t 

x p(^-f,-W^) (j|) j G (^-3 W+ , )(t)?))! 

\9-r2) r l VI) R \R^~ 1+R-t T 

(- s - £ -i,^-| )(r?) x G(g+ , +i7l+ ,_ |)(t ^ ); 



) 



x P 



(LI) 
(L2) 

(Jl) 

(J2) 



where R = f ^1 - 2-qt+t 2 . 

Next let us introduce the double generating function, that is the generating function of the 
generating functions Gi(t,rj; A): 



£=0 



For LI and L2 cases, the explicit forms are: 

(LI): G(s, t, rj; g) = 2 a ~i ((2 - t)y/l - t + ty/1 - t - 4s) 



1 5 



(L2) : d (td t ,g)g(s,t,r,;g) = [g+ \ - t ^-l+fi ri 



(1-t + s) 

which are obtained by using the two shifted generating functions (|E.9P and ()E.10|) in Ap- 
pendix EU It is a good challenge to derive the double generating functions for the Xi Jacobi 
polynomials. 



10 Three term recurrence relations 

Three term recurrence relations are one of the most fundamental characteristics of the ordinary 
orthogonal polynomials of one variable. Obviously the exceptional orthogonal polynomials do 
not satisfy these relations. Nevertheless, being deformations of ordinary orthogonal polynomials, 
the Xi polynomials are expected to retain certain reminiscent properties of the three term 
recurrence. 

Here we present a simple modification of the three term recurrence relations valid for the Xe 
polynomials. Its relevance is, however, as yet unclear. Let us denote the three term recurrence 
relation for the Laguerre or the Jacobi polynomials as 



rjP n ( m A) = A n (\)P n+1 ( V ; A) + B n (X)P n ( m A) + C n (X)P n ^ 1 (r ] ; A), 
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in which the explicit forms of the coefficients A n , B n and C n can be read from (lET5|) and (|ETl8|) 
in Appendix IE . 1 1 and Appendix lE.21 As a substitute of the above three term recurrence relations, 
we expect that a certain element in the degree £ + n + 1 invariant polynomial subspace Vi n +i 
(|6.3p . (|6.4p . which is related to r}P^ n {rj;\'), to be expressed in terms of Pi jn+ i(r}; A), P^^ n {ri;X) 
and P^ n _i(r/;A). From (|2.8p this can be achieved as 

E e> x[vPn(v, A')] = A n (X')d (n + 1, \)P^ n+l { m A) + B n (\')d (n, X)P e , n ( V ; A) 
+ C n (\')d (n - 1, A)P^ n _i(r?; A), 

where A' = A + M + <5. 

11 Zeros of polynomials 

The zeros of orthogonal polynomials have always attracted the interest of researchers. In the 
case of Xi polynomials P^ n (7/;A), it has n zeros in the domain where the weight function is 
defined, that is (0, oo) for the LI and L2 polynomials and (—1, 1) for the Jl and J2 polynomials. 
See for example, Section 5.4 of [19]. The behaviour of these zeros are the same as those of 
other ordinary orthogonal polynomials. This is guaranteed by the oscillation theorem of the 
one-dimensional quantum mechanics, since P£ n (j]; A) are obtained as the polynomial part of the 
eigenfunctions of a shape invariant quantum mechanical problem [HE]. 

Here we discuss the location of the extra £ zeros of the exceptional orthogonal polynomials, 
which lie in various different positions for the different types of polynomials. So far we have 
verified by direct calculation for lower £ and n: The £ extra zeros of LI polynomials are on 
the negative real line (— oo,0). Those of the L2 A^odd polynomials are 1 real negative zero 
which lies to the left of the remaining — 1) pairs of complex conjugate roots. The L2 A^ :even 
polynomials have \£ pairs of complex conjugate roots. For L2, these £ additional roots lie to 
the left of the n real zeros. 

The situations for the Jacobi polynomials are a bit more complicated. The Jl A^ :0 dd 
polynomials have 1 real negative root which lies to the left of the remaining \{£ — 1) pairs of 
complex conjugate roots with negative real parts. The Jl X£ :even polynomials have \£ pairs of 
complex conjugate roots with negative real parts. For Jl, some of the complex roots have real 
parts between —1 and 0. The J2 A"^ :0 dd polynomials have 1 real positive root which lies to the 
right of the remaining ^(£—1) pairs of complex conjugate roots with positive real parts. The J2 
A^ : even polynomials have \£ pairs of complex conjugate roots with positive real parts. For J2, 
some of the complex roots have real parts between and 1. 

12 Summary and comments 

We have given an in-depth study of the properties of the exceptional (Xi) polynomials discovered 
recently in [TJ [2J [IS]. Our main focus is the derivation of certain equivalent but much simpler 
looking forms of the Xg polynomials. The derivation is based on the analysis of the second order 
differential equations for the Xi polynomials within the framework of the Fuchsian differential 
equations in the entire complex plane. These new forms of the X% polynomials allow easy 
verification of the actions of the forward and backward shift operators on the Xi polynomials, 
and provide direct derivation of the Rodrigues formulas and the generating functions. The 
structure of the invariant polynomial subspaces under the Fuchsian differential operators is 
elucidated. The bases of the invariant polynomial subspaces provide a simple substitute of 
the three term recurrence relations. The Gram-Schmidt construction of the Xi polynomials 
starting from the above bases is demonstrated with the help of an integration formula. The 
properties of the extra zeros of the X^ polynomials are discussed. Some technical details are 
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relegated to the Appendices. The proof of the equivalence of the new and original forms of the 
Xi polynomials is given. Simple proofs of the forward and backward shift operations are shown. 
The integration formula is elementarily proven. Various fundamental formulas of the Laguerre 
and Jacobi polynomials are supplied for easy reference. 

Let us mention that the same method, deformation in terms of a degree £ eigenpolynomial, 
applied to the discrete quantum mechanical Hamiltonians for the Wilson and Askey-Wilson 
polynomials produced two sets of infinitely many shape invariant systems together with excep- 
tional (Xi) Wilson and Askey-Wilson polynomials (£ = 1,2, . . .) [IB]. It will be interesting to 
carry out the same analysis in these discrete cases. 

Finally, concerning the issue of global solutions of Fuchsian differential equations, we would 
like to make a comment on the well-known theorem by Heine-Stieltjes |20| . It asserts the 
existence of a polynomial solution for the differential equation 

d 2 d 
A(x)-^y(x) + B(x)—y(x) + C(x)y(x) = 0, 

in which only two coefficient functions A(x) and B(x) are specified. They are degree p + 1 
and p polynomials, respectively. In this case a polynomial C(x) of degree p — 1 is not given at 
the beginning but is determined so that the equation admits a degree n polynomial solution. 
Thus the problem setting is not a proper Fuchsian differential equation and the process of 
determining C(x) and the polynomial solution is purely algebraic. We consider these polynomial 
solutions do not qualify to be 'global solutions' of the ordinary Fuchsian equations. 

After this paper was arXived, some of the results were re-derived in terms of the Darboux- 
Crum transformations |17l [T8] . 



A Equivalence of different forms of polynomials 

The exceptional Jacobi polynomial for the trigonometric DPT presented in PQ, namely J2 ex- 
ceptional polynomial, is 

p t \\ dci (c( j.1 2n(-g + h + £-l)^ 1 ( V ;g,h + 2) 

A) = &(r); g + l,h + l) + 



(-g + h + 2£-2)(g + h + 2n + 21 - 1) 
n(2h + 4e-3)&- 2 (v,9 + l,h + 3y 



(2g + 2n + l)(-g + h + 2£- 2) 



P n ( V ;\ + £d) 



(-g + h + £- l)(2g + 2n + 2£- IV . . , /)s\ 

+ (2g + 2n + l)(g + h + 2n + 2£-l) ^ ff < H + + ^ 

In this paper we have presented it in a much simpler form in (|2.3p (J2) 

P i>n (n; A) = r ((g + g + l,h+ l)P n ( V ; g + £+l,h + £-l) 



n + g + 



(1 - v )b( m g, h)d v P n ( V ; g + £ + l,h + £-l)). 



In the following we write them as P° T ^(r];X) and P°° w (r/;A), respectively and show that they 
are in fact equal by using various identities of the Jacobi polynomials. For lower £ and n the 
equality can be verified by direct calculation. 

We fix £ and use new parameters a and (3 instead of g and h, 
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By using the forward shift relation for the Jacobi polynomial (|E.15p . the polynomials £e and P^ n 
are expressed as 

A) = Pt^'^Hv), A) = Plf- 2 \v\ 



n(P + £-l)p£f\ V ) 



{a + n- & + l){2l- a + (3 -2) 



(21 - a + (3 - 2)(2n + a + (3) 



(l_ a + /3 _l)( a + n ) ( _ a/?) (Qi/3) 
+ (a + n _^ + l )( 2n + a + /3)^- 1 (A ' lj 

- i(n + a + /3 + 1)(1 - ^P^-^^^ 2 ^)) • (A.2) 

Here we provide the proof for the equivalence of the two expressions (lA.ip and (|A.2p step by 
step: 

,or g , . x , « /(^ + m ~ v)P£? S (ri) + (a + 1)(1 + viPilT^Hv) 



Pin g (V,X) 



-21 



2n(£-a + f3-l)Pl«> p> ( V ) + 



(2^ - a + /3 - 2)(2n + a + /3) (a + n - ^ + 1)(2£ - a + - 2) 

2(1 - i)pj:r%) - (-a + g + 2£ - 2 )(i + t,)(p£?%) - pjrr 1,g+1) ^)) 

X -2(/3 + £-l) 

(a + n — £ + l)(2n + a + p) 
x 2nP^'%) - (a + g + 2n)(l - rj^fHjj) - P^'^jr,)) 

2(a + n) 



n + a + 1 

\ r \pt-i + i)\ < y + Ti)±'i_ x 



1 f , , ^ ^ + /?)(! - ^)^-T P; fa) + (a + 1)(1 + r?)P;_-r ™ W 

(a — I + 1)- 



a + n-< + ir ' -21 

x ((n + a + l)Pi Q '%) - (a + l^ 1 ^ 1 )^ 

( S - 

a + r; 

+ Pi~ a ~ X,P ~ X \ri){{n + a + lji**^ fa) - (a + 1)P^ +1 ^ 



_^( (a _, +1) p(--.« { , )P <„ +1 ,,- 1)w 



Properties of the Exceptional (Xp) Laguerre and Jacobi Polynomials 



19 



® 1 ((a - / + l)^— 2 ^ (^P^" 1 ) fo) 



a + n - £ + 1 

+ Pt^'^Hv^in + a + (3 + 1)(1 - ?? )P£f'%)) = PT(^ A )' 



where we have used (|R23l) . (lET27j) and (p^ijl in (i), (jET23|) and (lE^6|) in (ii), and flR25]) in (iii). 

The other exceptional polynomials, the Jl case is obtained from the above J 2 case by (|E.14|) . 
The equality of the exceptional Laguerre polynomials given in [TJ [2] and those given in this 
paper (|2.1f) will not be given here, since the LI and L2 cases are obtained from the Jl and J2 
cases by the limit (|E.28|) . 



B Forward and backward shift relations 

Here we provide proofs for the forward (]4.3p and backward (|4.4p shift relations which apply 
equally for the four types of Xj> polynomials. The method is elementary based on various 
identities of the Laguerre and Jacobi polynomials. 

The forward shift relation (14.3P is equivalent to a polynomial identity (|5.ip . By using do (n — 
1, A + 8) = do(n, A), / n (A + 6) = f n (X) and (|2.23p . it is easy to show that r.h.s. of f)5. 1 1) can be 
factorised d v P n (rj\ A + £8 + 8) x (•••). The forward shift relation is thus equivalent to (•••) = 0, 
namely, 

= di(A)&07; A + S) 2 - di(A + 6)&(ri; AJ&fa A + 28) - d v d 2 ( V )^( V ; A)^(r?; A + 8) 

+ d 2 (v)Uv, X)9r,tt(m A + 5) - d 2 (r,)d n Ce(v, X)Uv, A + 5). (B.l) 

This is a polynomial identity of degree 2^ and it is quadratic in the Laguerre/ Jacobi polynomials. 
This identity can be proven elementarily by using g?i(A) — d\(X + 8) = d v d 2 (ij) and (I2.20P to 
eliminate £^(77; A + 28) and ^^(ry; A). 

The backward shift relation (|4.4p is equivalent to a polynomial identity ()5.2p . By using 
c?o( n — 1, A + 5) = do( n ) A) and 6 n _i(A) = — 2n, the r.h.s. of (|5.2p becomes 

P n fa A') x (• • • ) + fyP^; A') x (• • • ) + P^fa; A' + 8) x (• • • ) 
+ a„P„_i(7 ? ; A' + 8) x (• • • ) + d 2 P n ^(r]; X' + 8) x (■■■), 

where A' = A + £8 + 8. The above expression can be reduced to P n _i(?]; A' + 8) x X(rj) + 
d v P n -i(r];X' + 8)xc 2 {v)y{v), by using ([233]) . CLU, (|2TT7p and the relations (ci(r/; X + 8 + 8) - 
ci(r]; X))d 2 (r]) = 2c 2 (if)d rj d 2 (if) with £ n (X) = / n (A)6 n _i(A). Up to an overall normalization, this 
Y{rj) is just the r.h.s. of (jB.ip . so it vanishes. Hence the backward shift relation is equivalent to 
X (rj) = 0, namely, 

= dfa, W&fa A + <5) 2 - ci(77, A + £S)d 1 (X + S)&( V ; AJfcfa; A + 25) 
- c 2 {rf)di{X + 5) (fcfa A)$,&fa; A + 28) - d^m \)Si(v, X + 25)) 

+ ^(A'^W&fa; V&im a + 5). (b.2) 

This is a polynomial identity of degree 2£ + 1 and it is quadratic in the Laguerre/ Jacobi poly- 
nomials. With the help of (|2.20p and (|2.18p . it is elementary to show that the r.h.s. of (|B.2p 
becomes £e(j]l A + 8)di(X)d 2 (rj) x (•••). Up to an overall normalization, this (• • • ) part is just 
(l.h.s. of (|2.2ip ) — (r.h.s. of (|2.2ip ). so it vanishes. This concludes the proof of the backward 
shift relation. 



20 



C.-L. Ho, S. Odake and R. Sasaki 



C Proof of the integration formula 

Here we present a simple proof for the integration formulas (|7.ip or (|T.2|) . We evaluate 
(BeMri)], Et,x[q(v)]) . x = ! {diWUv, A + S)p(r,) - d 2 (vMv, X)d v p(r,)) 



< (di(\)b( m A + 6)q(ri) - MvMm A)3,gfo)) W j%*M 6) dr ) 

+ 1 1 frfaA) J W^A-hM) I ^> ^ 

( = } y (di(A)d 3 (A + £)p(v)q(v)W(r,; A') + ^M^M^fa A' + 5))dr,, 

where A' = A + £8 + S. In (i) we have integrated by part and in (ii) we have used 
'di(A)6fo; A + d)\ 2 8n(W(m A + ^f^ +S U 2 ( V )) 



te(ril\) ) W(r l ;X + £S) 

= d 1 (X)d 3 (X + £d,£)^^-, 

which is shown by using (|2. 18j) and (|2.20p to eliminate ^(r?; A + S). 

D Hyperbolic DPT potential 

Here we provide a brief summary of the properties of the Xi Jacobi polynomials related to the 
deformed hyperbolic DPT potential. They are of the J2 type. The I = 1 case was introduced 
in [11] and the general I case was studied in pQ. In contrast to the radial oscillator and the 
trigonometric DPT potentials, the undeformed hyperbolic DPT potential allows only a finite 
number of square integrable polynomials P n (rj, A), n = 0, 1, . . . ,n B = [\{h — <?)]', where [a;]' 
denotes the greatest integer not equal or exceeding x. The situation is the same for the excep- 
tional polynomials, Pj> jn (r],\), n = 0, 1, . . . ,ng — I. Except for this point, the arguments for 
the hyperbolic DPT are the same as the radial oscillator and the trigonometric DPT potential 
cases. In particular, the much simpler looking new forms of the polynomials are also equivalent 
to the original forms of the polynomials given in pQ. We present various data: 

x^(g,h), h> g >o, <^ f (i,-i), p n (v;g,h) *p!T*-^\t,), 

tU{Wi9,h) = P e (77), n<n B , £<n B , 

Pi,n(v,X) (USD, n = 0,l,...,n B -£, 

do(n, A) = f n + g+ \, d= f (l,l), di{X) d = g+\, d 2 (r]) d = 1 - 77, 
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W(V,X) = 25^TT(^- ly-Hv + l)- h -3, 1< 77 < OO, 

, def r(n + g + ^)r(/i - g - n + 1) 
n ^ ' ~ 2n\(h-g-2n)T{h-n+\Y 

he,n{9,h) = h n (g + £,h - £)- 



(n + g+±)(h-n-l+± 



ci(r],\) d = g + h+(g- h+ l)r], c 2 (r/) = f rf - 1, 

61(77, X,£) = h-g-2£ + l-(g + h) V , 

S n (X) = An{h -g- n), E t {\) = A£{g + h + l-£), 

d 3 (X,£) = h-£+l, c T = A, f n {X) = 2{n + g-h), 6 n _i(A) d = -2n, 

n-l 



Uv,X) yrfd , a , _&fa;A + (fc + l)^ 

fa- l)^(r? + l)-^iV 

^(77; A + nS) 
d (td t , X)G e (t, n;\) = J2 t n d (n, X)P e , n ( V ; A) 

n 



&(7?;A + (n + l)<5) 



U/~l~ KV) r \r-t i + R-t "+" \ + R + t) r i VI) 

V/,„ = f Span[r/ fc ( 5 + ±)&fa; 5 + 1, & - 1) - *(1 - rfirf- 1 b( m g, h); k = 0, 1, . . . , n], 

Pn(v) = (V ~ !)" for 



us—I 

In the expression of the generating function, the summation range is extended to infinity, ~~ ^ 

n=0 

, formally. The forward and backward shift relations (|5.ip and (|5.2[) in terms of the Jacobi 

n=0 

polynomials are the same as the J2 case. 

E Summary: properties of the polynomials 

E.l Laguerre polynomials 

• Definition (expansion formula) 

= + k + l)„_ fc x fc . (E.l) 

n! fc! 

fc=o 

• Forward and backward shift relations 

d x Lt\x) = -L { n a _ + 1 1 \x), (E.2) 
xd x L^\x) + (a + 1 - xjL^Cx) = nL^(x). (E.3) 

• Differential equation 

x8%Lto (x) + (a + 1 - x^L^ (x) + raL&*> (x) = 0. (E.4) 
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• Three term recurrence relation 

xL^{x) = -(n + l)Ll%(x) + (2n + a + l)L^{x) - (n + (E.5) 

• Rodrigues formula 

1 1 / r/ \ n 

4") ( z ) = ( (e^x^) . (E.6) 

• Orthogonality (a > — 1) 

/ dxe- x x a Li a \x)L^(x) = -T(n + a + l)8 nm . (E.7) 
Jo n - 

• Generating functions 



tx 



G^(t,x) t f ^f4 a '(x) = (E.8) 



n=0 



G { «\t,x) d ^yt n L^ +n \x) = . a , (E.9) 



OO 



= f ^t n 4 Q " n )(x) = (1 +t) a e~ tx . (E.10) 



n=0 



Formulas (|E.9[) and (jE.lOp can be derived by using (|E.2[) and (|E.4|> . 

• Identities 

LM(x)-L£-V(x)=L ( ?l 1 (x), (E.ll) 

xL^x) - aL^x) = -nL^(x), (E.12) 

which can be verified elementarily based on fjE. 1 1) . 

E.2 Jacobi polynomials 

• Definition (expansion formula) 

• Parity 

pW'(-,) = (-l) n p(M(x). (E.14) 

• Forward and backward shift relations 

d x Pt P \x) = l(n + a + P + l)P^ +1 \x), (E.15) 
(1 - x 2 )d x P^ +1 \x) + (p- a -{a + P + 2)x)P ( n a X l ^ +l \x) = -2nP^(x). (E.16) 

• Differential equation 

(1 - x 2 )d 2 x P^\x) + (/?-a-(a + /3 + 2)x)d x P^\x) 

+ n{n + a + p + l)P^\x) = 0. (E.17) 



Properties of the Exceptional (Xe) Laguerre and Jacobi Polynomials 



23 



• Three term recurrence relation 

2( n + l)(n + a + f} + l)P^?\x) _ a 2 )P^\ x ) 



x 



(2n + a + /3 + l)(2n + a + /3 + 2) (2n + a + /3)(2n + a + /3 + 2) 

2(n + a)(n + /3)P^f ) (x) 
(2n + a + /3)(2n + a + /3 + 1) ' 1 ' ' 



Rodrigues formula 



p(^)(x) = L_^_ _i (JL) ((l-x) n+Q (l + x) n+/3 ). (E.19) 

n v ; 2"ra! (1 - x) a (l + x)P \dx J yK J v 1 } y 1 



• Orthogonality (a,/3 > — 1) 
i 



dx(l - x) a {l + xfP^\x)P^\x) 

_ 2"+/ j + 1 r(n + q + l)r(n + /3 + l) 
" n! (2n + a + J) + l)T(n + a + /3 + 1) 
Generating function 



(E.20) 



G^\t,x) ^TfpM(x) = — - j- ™, (E.21) 

v ' y ^ n w P(l + P-i) Q (l + P-K)^ v ; 

n=0 

R = y/l - 2xt + t 2 . 

• Identities 

2(n + /SjP^-^Cx) - 2/3P^- 1 ^(x) = (n + a + /3)(1 + s)p£f +1 \x), (E.22) 

(n + - x)p£+ 1 ^ ) (x) - q(1 + x)P^f +1) (x) = -2nP( a - l ^(x), (E.23) 

2(a + n)P^f \x) + (a + /3 + 2n)(l - x)(P^(x) - P^+^-^x)) = 2nP^(x),(E.24) 

2(n + a)P^" 1 ' /3 )(x) -2aP^' /3 - 1 )(x) = -{n + a + /3)(1 - x)pjff' P) (x), (E.25) 

(n + a)(l + x)P^f +1) (x) - 0(1 - x)P^+ 1,/3) (x) = 2nP^ /3 - 1 )(x), (E.26) 

-2( ( S+n)P^f ) (x) + (a+/3+2n)(l + x)(P^(x) - P^ a ~ l ^ +1 \x)) = 2nP^\x). (E.27) 

Equations (|E.22p - (|E.24p can be verified elementarily based on ()E.13p . By using (|E.14p . equa- 
tions (|E.25|> — (|E.27|) are obtained from equations (|E.22[> — (|E.24|) respectively. 

• Limit to the Laguerre polynomial 

mn P^(l - f ) = L^{±x). (E.28) 

Various formulas for the Jacobi polynomials reduce to those for the Laguerre polynomials in 
this limit. 



Acknowledgments 

This work is supported in part by the National Science Council (NSC) of the Republic of China 
under Grant NSC 96-21 12-M-032-007-MY3 (CLH), and in part by Grants-in-Aid for Scientific 
Research from the Ministry of Education, Culture, Sports, Science and Technology, No. 19540179 
(RS). Part of the work was done during CLH's visit to the Yukawa Institute for Theoretical 
Physics (YITP), Kyoto University, and he would like to thank the staff and members of YITP 
for the hospitality extended to him. RS wishes to thank National Taiwan University and National 
Chiao-Tung University for the hospitality extended to him during his visits in which a part of 
the work was done. 



24 



C.-L. Ho, S. Odake and R. Sasaki 



References 

Odake S., Sasaki R., Infinitely many shape invariant potentials and new orthogonal polynomials, 
\Phys. Lett. 5| 679 (2009), 414-417, larXiv: 0906. 01421 

Odake S., Sasaki R., Another set of infinitely many exceptional (Xe) Laguerre polynomials, \Phys. Lett.~B\ 
684 (2010), 173-176, larXiv:0911.3442l 

Infeld L., Hull T.E., The factorization method, \R~ev. Modern Phys~] 23 (1951), 21-68. 

Cooper F., Khare A., Sukhatme U., Supersymmetry and quantum mechanics, \Phys. Rep. | 251 (1995), 267- 



385, hep-th/9405029 



Darboux G., Sur une proposition relative aux equations lineaires, C.R. Acad. Paris 94 (1882), 1456-1459. 
Darboux G., Legons sur la theorie des surfaces, Vol. 2, 2nd ed., Gauthier-Villars, Paris, 1915, 210-215. 
Poschl G., Teller E., Bemerkungen zur Quantenmechanik des anharmonischen Oszillators, Z. Phys. 83 
(1933), 143-151. 

Bochner S., Uber Sturm-Liouvillesche Polvnomsvsteme, [Math~Zl 29 (1929), 730-736. 

Gomez-Ullate D., Kamran N., Milson R., An extension of Bochner's problem: exceptional invariant sub- 
spaces, \J~Approx. Theory\ l62 (2010), 987-1006. larXiv:0805. 33761 

Gomez-Ullate D., Kamran N., Milson R., An extended class of orthogonal polynomials defined by a Sturm- 
Liouville problem, [J Math. Anal. Appl.\ 359 (2009), 352-367. larXiv:0807.3939l 

Gomez-Ullate D., Kamran N., Milson R., Supersymmetry and algebraic Darboux transformations, 
| J. Phys. A: Math. Gen. 1 37 (2004), 10065-10078, |nlin.SI/0402052| 

Gomez-Ullate D., Kamran N., Milson R., Quasi-exact solvability and the direct approach to invariant sub- 
spaces, \J7Thys. A: Math. Gen. | 38 (2005), 2005-2019, [nlin.SI/0401030| 

Gendenshtein L.E., Derivation of exact spectra of the Schrodinger equation by means of supersymmetry, 
JETP Lett. 38 (1983), 356-359. 

Quesne C, Exceptional orthogonal polynomials, exactly solvable potentials and supersymmetry, 
\J. Phys. A: Math. Theor\ A\ (2008), 392001, 6 pages. larXiv:0807.4087l 

Bagchi B., Quesne C, Roychoudhury R., Isospectrality of conventional and new extended potentials, second- 
order supersymmetry and role of TT symmetry, \Pramana J. Phys.\ 73 (2009), 337-347, larXiv:0812. 14881 

Crum M.M., As sociated Sturm-Liouville systems, | Quart. J. Math. Oxford Ser. (2j\ 6 (1955), 121-127, 



lphysics/9908019] 

Quesne C, Solvable rational potentials and exceptional orthogonal polynomials in supersymmetric quantum 
mechanics, \SIGMA\ 5 (2009), 084, 24 pages, larXiv:0906.2331l 

Tanaka T., A/"-fold supersymmetry and quasi-solvability associated with JG-Laguerre polynomials, 
| J. Math. PTvys~\ 51 (2010), 032101, 20 pages. larXTv:0910.0328l 

Odake S., Sasaki R., Infinitely many shape invariant potentials and cubic identities of the Laguerre and 
Jacobi polynomials, | J. Math. Phys.\ 51 (2010), 053513, 9 pages, larXTv:0911.1585l 

Odake S., Sasaki R., Infinitely many shape invariant discrete quantum mechanical systems and new ex- 
ceptional orthogonal polynomials related to the Wilson and Askey- Wilson polynomials, \Phys. Lett. B 682 
(2009), 130-136. larXiv:0909.3668l 

Gomez-Ullate D., Kamran N., Milson R., Exceptional orthogonal polynomials and the Darboux txaxmiox- 
mation, \J~Phys. A: Math. Theor. || 43 (2010), 434016, 16 pages. larXTv:1002. 26661 

Sasaki R., Tsujimoto S., Zhedanov A., Exceptional Laguerre and Jacobi polynomials and the corresponding 
potentials through Darboux-Crum transformations, J. Phys. A: Math. Theor. 43 (2010), 315204, 20 pages, 
larXiv:1004.47TT1 

Andrews G.E., Askey R., Roy R., Special functions, Encyclopedia of Mathematics and its Applications, 
Vol. 71, Cambridge University Press, Cambridge, 1999. 

Heine H.E., Theorie der Kugelfunctionen und der verwandten Functionen, Berlin, 1878. 

Stieltjes T.J., Sur certains polynomes qui verifient une equation differentielle lineaire du second ordre et sur 

la theorie des fonctions de Lame, Acta Math. 6 (1885), 321-326. 



